February 1, 2008 21:18 WSPC/Trim Size: 9in x 6in for Proceedings Vill 



o 
o 

(N 



X 



STATISTICAL HADRONIZATION OF CHARM IN HEAVY 

ION COLLISIONS 



, A.P. KOSTYUK 

Institut fiir Theoretische Physik, Universitdt Frankfurt, Germany 

and 

^ ' Bogolyubov Institute for Theoretical Physics, Kyiv, Ukraine 

m ■ 
(N 

Production of open and hidden charm hadrons in heavy ion collisions is considered 
' within the statistical coalescence model (SCM). Charmed quark-antiquark pairs are 

, assumed to be created at the initial stage of the reaction in hard parton collisions, 

f"/^ 1 The number of these pairs is conserved during the evolution of the system. At 

f**^ ' hadronization, the charmed (anti)quarks are distributed among open and hidden 

charm hadrons in accordance with laws of statistical mechanics. Important special 
cases: a system with a small number of charmed quark-antiquark pairs and charm 
hadronization in a subsystem of the whole system are considered. The model 
^ • calculations are compared with the preliminary PHENIX data for J/psi production 

^ ' at RHIC. Possible influence of the in-nuclear modification of the parton distribution 

(-H ^ functions (shadowing) on the SCM results is studied. 



The thermal hadron gas (HG) model explains successfully the chem- 
' ical composition of light hadrons produced by colliding nuclei in a wide 

range center-of-mass energies^. A large body of experimental data can be 
well fitted with only three^ free parameters: the temperature T, baryonic 
chemical potential /i^ and volume V at the point of chemical freeze-out. 
This suggested the idea^ to see, whether this model can be applied also to 
charm particles. A straightforward generalization, when particles contain- 
ing heavy (anti)quarks are treated in the same way as light hadrons, could 
hardly be successful. The estimated relaxation time for charm production 
and annihilation in HG or even in a deconfined medium appeared to be 
much larger than the total lifetime of the thermal system. Therefore, it has 
been generally believed that the production mechanism of heavy-flavored 
hadrons is completely different from that of light ones. For instance, the 
standard picture of charmonium production in nucleus-nucleus collisions 



^Sometimes one more parameter, the strangeness suppression factor 7s is introduced. 



1 



February 1, 2008 21:18 WSPC/Trim Size: 9in x 6in for Proceedings 



Vill 



2 



assumes that charmonia are created exclusively at the initial stage of the 
reaction in primary nucleon-nucleon collisions. During the subsequent evo- 
lution of the system, the number of hidden charm mesons is reduced because 
of absorption of pre-rcsonance charmonium states in the nuclei (the nor- 
mal nuclear suppression) , interactions of charmonia with secondary hadrons 
(comovers), dissociation of cc bound states in the deconfined medium. The 
last mechanism was expected to be cspcicially strong and charmonia were 
proposed to be used as a probe of the state of matter created at the early 
stage of the collision. It was found that the J/ ^ suppression with respect 
to Drell-Yan muon pairs measured in proton-nucleus and nucleus-nucleus 
collisions with light projectiles can be explained by the normal nuclear 
suppression alone^. In contrast, the NA50 experiment with a heavy pro- 
jectile and target (lead-lead) revealed essentially stronger Jf-ip suppression 
for central collisions^. This anomalous J /ijj suppression was attributed^ to 
formation of quark-gluon plasma (QGP). In the same time, purely hadronic 
scenarios still cannot be excluded''. Despite of quite successful agreement 
with the J/tp data, the standard scenario seems to be in trouble explaining 
the ijj' yield. The recent lattice simulations* suggest that the tempera- 
ture of if)' dissociation Tci{'4>') lies far below the dcconfinemcnt point^ Tc- 
Td{tp') ~ 0.1-0.2rc. Therefore, not only the quark-gluon plasma, but also 
a hadronic co-mover medium should completely eliminate i/j' charmonia in 
central Pb-|-Pb collisions at SPS. However, the experiment revealed a siz- 
able tp' yield (see, for instance^"). It was observed^-^ that ip' to J/tp ratio 
decreases with centrality only in peripheral lead- lead collisions, but satu- 
rates at sufficiently large nmnber of participants Np > 100. The value of 
the ratio in (semi)central collisions is approximately constant and equal to 
the ratio of the densities of these charmonium states in an equilibrium HG. 
Hence, on one hand, the production of charm cannot be thermal, because 
of the large relaxation time, on the other hand, the multiplicity ratio of 
different charmonium states is thermal. This paradox is resolved in the 
statistical coalescence model (SCM)^^. In this model, the charm quarks c 
and antiquarks c are created at the initial stage of A-l-A reaction in the 
hard parton collisions. This is similar to the standard approach. But like 
in the thermal modeP, the formation of observed hadrons with open and 
hidden charm takes place later at the hadronization stage near the point of 
chemical freeze-out. Production and annihilation of charm quark- antiquark 
pairs at all stages after the initial are neglected. Therefore the number c 
and c in the system may be very far from chemical equilibrium, but their 
distribution over different species of hadrons with open and hidden charm 
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is controlled by the lows of statistical mechanics. 

In the 'pure' SCM, which will be a subject of the present discussion, it 
is assumed that the hot strongly-interacting medium destroys all initially 
formed ('primordial') charmonia or even prevents their formation. A com- 
bined scenario, assuming that some of primordial charmonia do survive 
suppression in the medium, but additional hidden charm meson can be as 
well formed at hadronization, has been also considered^^. 

Within the grand canonical approach, the deviation of the amount of 
(anti)charm in the system from its equilibrium value can be taken into ac- 
count simply by multiplying the multiplicities of all open single (anti) charm 
species in equilibrium hadron gas by some factor 

{X) = jcNx{T, V) = 7c exp (cx^) Nx{T, ^lb, V). (1) 

Here Cx = ±1 is the charm of the spcc;ics X . The charm chemical potential 
fic is needed to keep the total balance of open charm and anticharm in the 
thermal system. At zero fic but nonzero baryonic chemical potential /ib, 
the total number of open charm in the equilibrium HG 

A^i = $^iVx(T,/x6,y) (2) 

X 

is not equal to that of anticharm 

Nx = J2^x{T,l^b,V). (3) 

X 

This is because, for instance, a positive value of /ih enhances the number of 
open charm baryons, that contribute to (2), and suppresses the number of 
antibaryons in (3). (The sums in (2) and (3) run over all species of single 
open charm and anticharm hadrons, respectively.) The value of fXc should 
be chosen to restore the balance: 

exp(/Xe/r)7Vi = exp(-Mc/T)7Vi. (4) 

A hidden charm meson contains a charm quark as well as an antiquark. 
Therefore, their thermal numbers are not influenced by the charm chemical 
potential. To take into account the deviation from the chemical equilibrium, 
one should multiply the thermal numbers of charmonia by 7^: 

(y)=7c'A^F(T,M6,F). (5) 

Provided that the average number of charm quarks and antiquarks 
(C) = (C) in the system is known, the factor 7^ can be found from the 
equation: 

{C) = 7c exp{fic/T)Ni + iINh • (6) 
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Here Nh is the total number of hidden charm particles in the system at 
chemical equilibrium; 

NH = ^NY{T,^i,,V). (7) 

Y 

(The sum runs over all charmonium states.) The thermal number of hidden 
charm is much smaller than that of open charm. Therefore, if the number of 
cc pairs in the system is not extremely large (does not exceed the chemical 
equilibrium value by a factor of several hundreds), the hidden charm term 
in (6) can be neglected and 7c is simply given by 

7c = ^^. (8) 

exp{pjT)N, 

The grand canonical approach does not exactly correspond to the phys- 
ical reality. It keeps balance only between the average quantities of charm 
and anticharm, while the partition function includes also configurations 
with nonequal numbers of charm quarks and antiquarks in the system. 
This does not influence the result, if the average number of cc pairs in the 
system is large {C),{C) » 1 . In the opposite case^^, however, it should 
be taken into account that strong interactions conserve charm (weak inter- 
actions can be neglected), therefore the number C of charm quarks in the 
system is always exactly equal to the number C of antiquarks. 

Let us first consider the situation, when the numbers of charm quarks 
and antiquarks are fixed (not necessary equal to each other). In this case, 
the partition function for particles with charm can be written as 

ZciT, M., V) = + (C^ (C3I)T + • • • 

The dots stand for terms corresponding to more than one hidden charm 

meson in the system and also for terms including double and triple charm 
baryons. These terms, as it was already mentioned above, can be safely 
neglected, unless C and C are extremely large. 

As far as Nh ^ NiNi (which is true in all practically interesting sit- 
iiations), the numbers of open charm and anticharm hadrons are approx- 
imately equal to C and C, respectively, i.e. almost all charm quarks and 
antiquarks hadronize into open charm particles. Only a very small frac- 
tion form charmonia. The total number of hidden charm can be easily 
calculated: 

^H)cc = Nh^^^^'^CC^ (10) 
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In reality, however, the number K of cc pairs is not fixed. It fluctuates 
from one event to another. The pairs are produced in independent collisions 
of nucleons whose number is large and the probability to produce a cc pair 
in a single collision is small. Therefore, the fluctuations should follow the 
Poisson law: 

Pi^C = C = K) = exp(-(cc-)^B(M)^^^^^^^, (11) 



where (cc)ab(6) stands for average number of cc pairs produced by nuclei 
A and B colliding at impact parameter h. 

The number of hidden charm averaged over all events at fixed centrality 
is given by the convolution of Eq.(lO) with the probability (11). The result 
reads: 

{H)AB{b) = {cc) AB(b){{cc) AB{h) + l)Tr-cr " (^^^ 

If one interested in the multiplicity of a particular charmonium species, 
one should replace the Nh in the numerator of (12) by the thermal multi- 
plicity of this species. It must be also taken into account that an additional 
contribution to multiplicities of low-lying charmonium states comes from 
decays of excited ones. 

Note that the ratios of the numbers of different charmonium states in 
SCM are exactly the same as in equilibrium HG: 

W)AB{b) ^ .^3. 

It agrees with the behavior of ip' to J /ip ratio in central Pb+Pb collisions 
at CERN SPS. This is not the case for peripheral events [Np < 100). It can 
be explained in the following way: SCM is valid if the momenta of charm 
(anti)quarks (not their number!) are thermalized. This is possible only in 
a sufficiently large system. 

The equation (12) gives the total (47r) multiplicity of hidden charm. In 
real experimental situation, however, measurements are made in a limited 
rapidity window. In the most simple case, when the fraction of charmonia 
that fall into the relevant rapidity window does not depend on the centrality, 
one can merely use Eq.(12) multiplied by some factor ^ < 1. This situation 
is likely to be relevant to charmonium production at CERN SPS, where the 
multiplicity of light hadrons, that determine the freeze-out volume of the 
system, are approximately proportional to the number of nucleon partici- 
pants Np at all rapidities. Indeed, SCM fits well the measured^ '^'^^ J/tjj 
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to DrcU-Yan ratio at CERN SPS (see Fig. 1, the details of the fit can be 
found in Ref.^^). 




Pb-Pb 1996 
Pb-Pb 1996 min. bias 
Pb-Pb 1998 min. bias 
Pb-Pb 2000 (set A) 
Pb-Pb 2000 (set B) 
Pb-Pb 2000 (set C) 
SCM 
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Figure 1. The dependence of the J/i/' to Drell-Yan ratio on the transverse energy at 
SPS. The vertical line shows the boundary of the applicability domain of the statistical 
coalescence model (SCM). 



At BNL RHIC, the situation is different: the total (47r) multipHcity of 
hght hadrons are still approximately proportional to the number of partic- 
ipants, while at midrapidity, it grows faster^''. The centrality dependence 
of charmonium production at different rapidities should in this case be also 
different. The formula (12) should be generalized, in such a way that the 
charmonium multiplicity not only in the whole system, but in a part of it, 
a subsystem (like a limited rapidity interval), could be calculated. 

Let, if c (or c) is present in the system, the probability to find it in the 
subsystem is ^ < 1. Then, if the number of cc pairs in the whole system 
is K, the probability to find C charm quarks in the subsystem is given by 
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the binomial law: 

We assume that the distributions of quarks and antiquarks are uncorre- 
lated,^ and the probability distribution of the number of antiquarks C is 
given by the same binomial law. 

If the number of pairs K is distributed in accordance with the Poisson 
law (11), the average multiphcity of charmonia in the subsystem is 

oo K K ^ 

wi'iw - E p^K) E ^(^1^) E ^(^1^) (15) 

= {cc)AB(b) {{cc)AB(b) + 1) 

IS \i \ 

Here the meaning of (cc)ab(6) is the same as in Eq.(12): the average to- 
tal (47r) number of cc pairs in an event with given centrality. In contrast, 
Nh, Ni and are related to the subsystem. These are the thermal mul- 
tiplicities of charmonia, open charm and anticharm calculated with the 
thermodynamic parameters fitted within the HG model to the chemical 
composition and multiplicity of light hadrons in the rapidity window under 
interest. 

In Fig. 2, the SCM result is compared to the preliminary PHENIX data^^ 
for J/ip production in the central rapidity interval at the top RHIC energy 
(^/s = 200 GeV). The values of freeze-out parameters were taken from 
Ref.^i; T = 177 MeV and /is = 29 McV. The dependence of volume on 
the centrality was chosen to reproduce the measured multiplicity of charged 
hadrons.^'' The value of {cc)AB(b) was calculated in the Glauber approach. 
The charm production cross section in nucleon-nucleon collisions was fixed 
at^" a^c^ = 650 /zb, while ^ was considered as a free parameter. The 
minimum value of x^/do/ — 0.75 is reached at ^ « 0.18. 

As is seen, the SCM dependence of the J/ip multiplicity at midrapid- 
ity per binary collision on the centrality is almost fiat for Np > 100 in 
contrast to the total J ftp yield , which is expected to grow in the same 
centrality region.^^ This is because the multiplicity of light hadrons^^ and, 
consequently, the hadronization volume at midrapidity grows faster with the 



''At this point, our approach differs from that of Ref.'^®, where a strong correlation C = C 
between the numbers of c-quarks and c-antiquarks is assumed. 
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Figure 2. The J/ip multiplicity per binary collision in the central unit rapidity interval 
at y/s = 200 GeV muliplyed by the J/V' — » e+e~ branching ratio versus the number of 
nucleon participants. 



centrality than the total vohimc. If the influence of the in-nuclear modifica- 
tion of the parton distribution functions (shadowing) on charm production 
at RHIC is essential, a decrease should be observed instead of saturation 
(see Fig.2). 

The centrality dependence of J/tp multiplicity in the QGP suppression 
scenario^'^ is also shown in Fig.2. The parameters of the model were fixed 
from SPS data^'^'^^ and then extrapollated to RHIC. As is seen, the ob- 
served yield can hardly be explained by survived primordial charmonia only. 
It seems that cc coalescence at late stages of the reaction is here a dominant 
if not the only charmonium production mechanism. 

In conclusion, the statistical coalescence model is consistent with the 
data on J/tp production in heavy nucleus collisions at SPS and RHIC ener- 
gies, while the standard scenario of J/tp dissociation in QCP predicts too 
strong suppression at RHIC that is not favored by the experimental data. 
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